Extension of the Conley-Zehnder Index and 
Calculation of the Maslov-Type Index Intervening 
in Gutzwiller's Trace Formula 



Maurice de Gosson 
Universitat Potsdam, Inst. f. Mathematik 
Am Neuen Palais 10, D-14415 Potsdam 
Current address: 
Universidade de Sao Paulo, Dep. de Matematica 
CEP 05508-900 Sao Paulo 
E-mail: maurice.degosson@gmail . com 

Serge de Gosson 
Vaxjo University, MSI, 
Vaxjo, SE-351 95 
E-mail: sergedegosson@gmail.com 

February 7, 2008 
Abstract 

The aim of this paper is to define and study a non-trivial exten- 
sion of the Conley-Zehnder index. We make use for this purpose of 
the global Arnol'd-Leray-Maslov index constructed by the authors in 
previous work. This extension allows us to obtain explicit formulae 
for the calculation of the Gutzwiller-Maslov index of a Hamiltonian 
periodic orbit. We show that this index is related to the usual Maslov 
index via Morse's index of concavity of a periodic Hamiltonian orbit. 
In addition we prove a formula allowing to calculate the index of a 
repeated orbit. 



Keywords: Gutzwiller trace formula, Conley-Zehnder index, Maslov in- 
dex. 

AMS Classification (2000): 37J05, 53D12, 53Z05, 81S30 



1 Introduction 



This paper is devoted to a redefinition and extension of the Conley-Zehnder 
index [3, 13] from the theory of periodic Hamiltonian orbits. This will allow 
us to prove two new (and useful) formulae. The first of these formulae 
expresses the index of the product of two symplectic paths in terms of the 
symplectic Cayley transform of a symplectic matrix; the second shows that 
the difference between the Conley-Zehnder index and the usual Maslov index 
is (up to the sign) Morse 's index of concavity. 

Besides the intrinsic interest of our constructions (they allow practical 
calculations of the Conley-Zehnder index even in degenerate cases) they 
probably settle once for all some well-known problems from Gutzwiller's 
[12] quantization scheme of classically chaotic systems. Let us briefly recall 
the idea. Let if be a Hamiltonian operator with discrete spectrum Ei,E2, ■■■■ 
One wants to find an asymptotic expression, for — > 0, of the "level density" 



oo ^ 



d{E) = ^5{E - Ej) = lim Tr G{x, x', E + ie) 

k=i ^ 

{G the Green function determined by H). Writing d{E) = d{E) + d{E) 
where d{E) is the Thomas-Fermi or smoothed density of states, and d{E) 
the "oscillating term", Gutzwiller's formula says that in the limit /I — > we 
have _ _ 

d{E)=dGntz{E)+0{h) 

where 

dGntz{E) = ^ Re > , = exp - * pdx . (1) 

, ^\detiS,-I)\ V^A J 

The sum in the right-hand side of the formula above is taken over all periodic 
orbits 7 with period of the classical Hamiltonian H, including their repe- 
titions; the exponent is an integer and and S^^ is the stability matrix of 7. 
Gutzwiller's theory is far from being fully understood (there are problems 
due to possible divergences of the series, insufficient error estimates, etc.). 
We will not discuss these delicate problems here; what we will do is instead 
to focus on the integers yU^, to which much literature has been devoted (see 
for instance [1, 4, 21] and the references therein). As is well-known, /i^ is 
not the usual Maslov index familiar from EBK quantization of Lagrangian 
manifolds, but rather (up to the sign) the Conley-Zehnder index of the orbit 
7: 

"Gutzwiller-Maslov index" = —(Conley-Zehnder index) 
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(see e.g. Muratore-Ginnaneschi in [18]; in [11] one of us has given a different 
proof using an ingenious and promising derivation of Gutzwiller's formula 
by Mehlig and Wilkinson [16] based on the Wcyl representation of meta- 
plectic operators). We notice that Sugita [23] has proposed a scheme for the 
calculation of what is essentially the Conley-Zehnder index using a Feyn- 
man path-integral formalism; his constructions are however mathematically 
illegitimate and seem difficult to justify without using deep results from 
functional analysis (but the effort might perhaps be worthwhile). 

This paper is structured as follows: 

1. Redefine the Conley-Zehnder index in terms of globally defined indices 
(the Wall-Kashiwara index and the Arnol'd-Leray-Maslov index); we 
will thus obtain a non-trivial extension v of that index which is explic- 
itly computable for all paths, even in the case det(>S' — /) = 0; this is 
useful in problems where degeneracies arise (for instance the isotropic 
harmonic oscillator, see [20] and the example of last section in the 
present paper); 

2. We will prove a simple formula for the Conley-Zehnder index of the 
product of two symplectic paths. We will see that in particular the 
index of an orbit which is repeated r times is 

^7r = f^i + \{'''- 1) sign Ms 

where 

Ms = \j{s+i){s-ir^ 

is the symplectic Cayley transform of the monodromy matrix S] 

3. We will finally prove that the Conley-Zehnder index v^^ of a non- 
degenerate periodic orbit 7 is related to its Maslov index m-y by the 
simple formula 

= rrij — Inert W-y 

where Inert W-y is Morse's "index of concavity" [17, 18], defined in 
terms of the generating function of the monodromy matrix. 

We close this article by performing explicit calculations in the case of the 
two-dimensional anisotropic harmonic oscillator; this allows us to recover a 
formula obtained by non-rigorous methods in the literature. 
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Notations 

We will denote by a the standard symplectic form on M^" = x 

<^{z,z') = {p,x') - {p',x) if z = {x,p), z' = {x'p') 
that is, in matrix form 

cf{z,z') = (^Jz,z') , J = 

The real symplectic group Sp(n) consists of all linear automorphisms S of 
M^" such that a{Sz, Sz') = a{z, z') for ah z, z' . Equivalently: 

S G Sp(n) ^ S'^JS = SJS^ = J. 

Sp(n) is a connected Lie group and 7ri[Sp(n)] = (Z, +). Wc denote by 
Lag(n) the Lagrangian Grassmannian of (R'^", a), that is: i € Lag(n) if and 
only ^ is a n-plane in M^" on which a vanishes identically. We will write 
£x = M^xO and^p = OxM^ (the "horizontal" and "vertical" polarizations). 

2 Prerequisites 

In this section we review previous results [7, 8, 10] on Lagrangian and sym- 
plectic Maslov indices generalizing those of Leray [14]. An excellent com- 
parative study of the indices used here with other indices appearing in the 
literature can be found in Cappell et al. [2]. 

In what follows {E,io) is a finite-dimensional symplectic space, dim£^ = 
2n, and Sp(£^, oj), Lag{E, oj) the associated symplectic group and Lagrangian 
Grassmannian. 

2.1 The Wall— Kashiwara index 

Let be a triple of elements of Lag{E,uj); by definition [2, 15, 24] 

the Wall-Kashiwara index (or: signature) T{i, £',£") is the signature of the 
quadratic form 

Q{z, z', z") = a{z, z') + a(z', /') + g{z\ z') 
on £ I" . The index r is antisymmetric: 

r(£,/,£") = -r(/,£,r) = -r(£,£",/) = -t(/',/,^); 



I 
-I 
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it is a symplectic invariant: 

T{Si, Si', St) = r(£, f ) for S G Sp(n) 
and it has the following essential cocycle property: 

r(£, - T{e, e", e") + r(^', i", e") - t(/, i", i'") = o. (2) 

Moreover its values modulo 2 are given by the formula: 

t(^,/,/') = " + dim£n/ + dim£'n£" + dim£"n£ mod 2. (3) 

Let {E, oj) = {E' e E", u' u"); identifying Lag{E', uj') Lag(£;", u") with 
a subset of Lag(£?, a;) we have the following additivity formula: 

T{h ^2 , i'l e 4 , 4' © 40 = Ti {h , e'l , O + t2 (£2 , 4 , 4) 

where ti and T2 are the signatures on Lag(£'',a;') and 'Lag{E" ,io"). 

The following Lemma will be helpful in our study of the Conley-Zehnder 
index: 

Lemma 1 (i) Ifini" = then T{i, £',£") is the signature of the quadratic 
form 

Q'{z') = uj(Pvee" z', z') = u{z', Pvt'e z') 

on m , where Pr^^" is the "projection onto I along i" and Pr^//^ = I — Vyu" 
is the projection on t" along £. (ii) Let {I, £',(.") he a triple of Lagrangian 
planes such that an£ = £r\e' + en £" . Then t{£, £', £") = 0. 

(See e.g. [15] for a proof). 

The index of inertia of the triple {£,£',£") is defined by 

Inert(^,/,/') = ^(r(£, f) + n + dim£n£' - dim^' n^" + dim^ n£); (4) 

in view of (3) it is an integer. When the Lagrangian planes £, £' , £" are pair- 
wise transverse it follows from the first part of Lemma 1 that Inert(£, £") 
coincides with the index of inertia defined by Leray [14] : see [7, 8] . 
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2.2 The ALM index 

Recall [7, 8] (also [10] for a review and [9] for calculations in the case n = 1) 
that the ALM (=Arnord-Leray-Maslov) index on the universal covering 
Lag^{E,Lo) of LsLg{E,io) is the unique mapping 

having the two following properties: 

• /X is locally constant on each set {(^oo ^to) ■ dimlCif = k} {0 < k < n); 

• For all ^ooj ^'oc5 ^oc ™ Lagoo(£^,w) with projections I, I" we have 

m(^oo,C) - m(^oo,C) + m(C,0 = r {£,£',£") (5) 
where r is the Wall-Kashiwara index on Lag(£^,a;). 
The ALM index has in addition the following properties: 

/"(^ocO = ^ + dim£n£' mod 2 (6) 

(n = ^ dim E) and 

Kf^'eocP'-'O = ^(^00,0 + 2(r - r') (7) 

for all integers r and r'; here /9 denotes the generator of Tri[Lag{E , uj)] = 
(Z, +) whose image in Z is +1. From the dimensional additivity property of 
the signature r immediately follows that if li^^o © ^2,00 and £[ 00 © ^2 00 are 
in 

Lag^(£;',a;') ©Lag^(^",^") C Lag^l^^,^) 

then 

M(^1,oo © £2,00, £'l,oc © ^2,00) = At' (^1,00, /l,oo) + M"(^2,oo, 4,00) (8) 

where and /j," are the ALM indices on Lag^(£",a;') and hag^{E" ,lo"), 
respectively. 

When {E, u) is the standard symplectic space (R^^, cr) the "Souriau map- 
ping" [22] identifies Lag(£^, a;) = Lag(n) with the set 

W(n, C) = {we U(n, C) : w = w'^} 



6 



of symmetric unitary matrices. This is done by associating to to ^ = n^p 
(n G U(n, C)) the matrix w = uu^; the Maslov bundle Lag^{n) is then 
identified with 

Woo{n,C) = {{w,e) ■.weW{n,C), detu; = e*^}; 

the projection tt^^^ : i^o i — > (■ becoming {w, 6) i — > w. The ALM index is 
then calculated as follows: 

• If £ n f = then 

M^oo,0 = - [e-e' + iTii.og{-w{w')-''] (9) 

TT 

(the transversality condition is equivalent to ~w{w')~^ having 

no negative eigenvalue); 

• If ^ n f 7^ one chooses any t" such that i r\ i" = I' r\ I" = and 
one then calculates fi{iao,i'^) using the formula (5), the values of 
m(-^oo)^^) m(-^to'^oo) being given by (9). (The cocycle property 
(2) of r guarantees that the result does not depend on the choice of 

see [7, 8]). 

2.3 The relative Maslov indices on Sp{E,io) 

We begin by recalling the definition of the Maslov index for loops in Sp(n). 
Let 7 be a continuous mapping [0, 1] — > Sp(n) such that 7(0) = 7(1), and 
set 7(t) = Sf. Then Ut = {StSt)~^^'^St is the orthogonal part in the polar 
decomposition of Sf. 

Ut e Sp(n) nO(2n,M). 
Let us denote by ut the image t{Ut) of Ut in U(n, C): 

LiUt) =A + iB if U= ^ ~f 

and define p{St) = det ut. The Maslov index of 7 is by definition the degree 
of the loop 1 1 — > p{St) in S^: 

m(7) = deg[t 1 — > det{i{Ut))] , < t < 1. 
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Let a be the generator of 7ri[Sp(£^, a;)] = (Z, +) whose image in Z is +1; if 
7 is homotopic to then 

m(7) = m(a'') = 2r. (10) 

The definition of the Maslov index can be extended to arbitrary paths in 
Sp(£', a;) using the properties of the ALM index. This is done as follows: let 
i = TT^^^iioo) e Lag{E,uj); we define the Maslov index of Soo G Sp^{E,Lo) 
relative to i by 

fJ'iiSoo) = KSood-ooJoo); (11) 

one shows (sec [7, 8]) that the right-hand side only depends on the projection 
£ of £oo: justifying the notation. 

Here are three fundamental properties of the relative Maslov index; we 
will need all of them to study the Conley-Zehnder index: 

• Product: For all ^oo, S'^, in Sp^{E,uj) wc have 

MSooS'^) = Ms^) + f^e{s'^) + Tie,si,ss'ey, (12) 

• Action of 7ri[Sp(n)]: We have 

l^eioTSoo) = MSoo) + 4r (13) 

for all r G Z; 

• Topological property: The mapping {Soo,i) i — ^ /^^(-S'oo) is locally con- 
stant on each of the sets 

{{Soo,i) :dimS£ne = k} C Sp^{E,uj) x Lag{E,uj) (14) 

{0<k< n). 

The two first properties readily follow from, respectively, (5) and (7). 
The third follows from the fact that the ALM index is locally constant on 
the sets {(^oo,^'oo) ■ dim£n^' = k}. Note that (13) implies that 

hence the restriction of any of the m to loops 7 in Sp(£^,ti;) is twice the 
Maslov index m{'y) defined above; it is therefore sometimes advantageous to 
use the variant of m defined by 

MSoo) = \{l^i{Soo) + n + dim(5^ n £)) (15) 
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where n = ^dim^^. We will call m£{Soo) the reduced (relative) Maslov 
index. In view of property (6) it is an integer; the properties of are 
obtained, mutatis mutandis, from those of /x^; for instance property (12) 
becomes 

meiS^S'^) = m^(-Soo) + me{S'^) + Inert(A S£, SS'£) 

where Inert is the index of inertia defined by (4). 

It follows from the cocycle property of the signature r that the Maslov 
indices corresponding to two choices £ and I' are related by the formula 

MSoc) - W'(5oo) = r{SiJ,f) - TiS£,Se',£'y, (16) 

similarly 

^^(-500) - meiSoo) = lnert{S£, £, £') - Inert(5^, S£', £'). (17) 

Assume that {E,u) = {E' E",u' u") and £' G Lag(£;',a;'), £" e 
Lag{E" ,uj"); the additivity property (8) of the ALM index implies that if 
S'^ G Sp^{E',u'), S':, G Sp^{E",u") then 

l^e'ee"{S'oo ® S'^,) = iJi,{S'^) + tJie,{S':^) (18) 

where ^^^{E' ,u)') © S\)^{E'\uj") is identified in the obvious way with a 
subgroup of S\y^{E,u); a similar property holds for the reduced index m^. 

3 The index v on Spoo(?T^) 

In this section we define and study a function u : Sp^{n) — > Z extending 
the Conley-Zehnder index [3] . We begin by recalling the definition and main 
properties of the latter. 

3.1 Review of the Conley— Zehnder index 

Let S be a continuous path [0, 1] — > Sp(n) such that S(0) = / and 
dct(S(l) — /) 7^ 0. Loosely speaking, the Conley-Zehnder index [3] counts 
algebraically the number of points in the open interval ]0, 1[ for which S(t) 
has 1 as an eigenvalue. To give a more precise definition we need some 
notations. Let us define three subsets of Sp(n) by 

Spo(n) = {S : det{S-I) = 0} 
Sp+(n) = {S : det{S-I) > 0} 
Sp-(n) = {5 : det(5-/) < 0}. 
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These sets partition Sp(n), and Sp''"(n) and Sp (n) are moreover arcwise 
connected; the symplectic matrices S'^ = —I and 



L 
L-i 



L = diag[2,-1,...,-1] 



belong to Sp"^ (n) and Sp~ (n) , respectively. 

Let us now denote by C±{2n, M) the space of all paths S : [0, 1] — > Sp(n) 
with S(0) = / and S(l) € Sp^(n). Any such path can be extended into 
a path E : [0,2] — > Sp(n) such that G Sp=^(n) for 1 < t < 2 and 
S(2) = S+ or S(2) = S~. Let p be the mapping Sp(n) — > S^, p{St) = 
det ut, used in the definition of the Maslov index for symplectic loops. The 
Conley-Zehnder index of E is, by definition, the winding number of the loop 
(poS)2 in S^: 

zcz(S) = deg[i ^ {pm))f, 0<t<2]. 

It turns out that icz(S) is invariant under homotopy as long as the end- 
point S = remains in Sp^(n); in particular it does not change under 
homotopies with fixed endpoints so we may view icz as defined on the subset 



Spl,{n) = {S^:dei{S-I)^0} 

of the universal covering group Sp(^(n). With this convention one proves [13] 
that the Conlcy-Zchndcr index is the unique mapping icz '■ Sp^(n) — > Z 
having the following properties: 



(CZi) Antisymmetry: For every Soo we have 

^cz('S'~^) = — jcz('S'oo) 

where is the homotopy class of the path t i — > S^^; 

(CZ2) Continuity: Let S be a symplectic path representing and T,' a 
path joining S to an element S' belonging to the same component 
Sp^(n) as S. Let S'^ be the homotopy class of S * S'. We have 

^cz('S'oo) = icz{S'^); 



(CZ3) Action of 7ri[Sp(n)]: 

icz{a''Soo) = icz{Soo) + 2r 

for every r G Z. 
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We observe that these three properties are characteristic of the Conley- 
Zehnder index in the sense that any other function zj-,^ : Sp^(n) — > Z 
satisfying then must be identical to icz- Set in fact 6 = iqz — i'cz- 
view of (CZ3) we have 6{a'^Soo) = (5(>S'oo) for all r £ Z hence 6 is defined 
on Sp*(n) = Sp+(n) U Sp-(n) so that d{Soo) = S{S) where S = Si, the 
endpoint of the path t 1 — > St- Property (CZ2) implies that this function 
Sp*(n) — > TL is constant on both Sp^(n) and Sp~(n). We next observe that 
since det S" = 1 we have det(S'~^ — 1) = det(S' — /) so that S and S~^ always 
belong to the same set Sp''"(n) or Sp^(n) if det(S' — /) 7^ 0. Property (CZi) 
then implies that b must be zero on both Sp'*"(n) or Sp~(n). 

Two other noteworthy properties of the Conley-Zehnder are: 

(CZ4) Normalization: Let Ji be the standard symplectic matrix in Sp{l). 
If is the path t — ^ e^*^'' (0 < t < 1) joining / to -I in Sp(l) then 
^cz,i (-51,00) = 1 (^cz,i the Conley-Zehnder index on Sp(l)); 

(CZ5) Dimensional additivity: if S'1,00 S Sp^(ni), <S'2,oo S Sp^(n2), ni + 
n2 = n, then 

^Cz('S'i,oo ffi S2,oo) = ^CZ,1 ('S'1,00) + ^CZ,2('S'2,oo) 

where iczj is the Conley-Zehnder index on Sp(nj), j = 1,2. 
3.2 Symplectic Cayley transform 

Our extension of the index icz requires a notion of Cayley transform for 
symplectic matrices. If S G Sp(n), det(S' — /) 7^ 0, we call the matrix 

Ms = y{S + I){S-I)-' (19) 
the "symplectic Cayley transform of S"' . Equivalently: 

Ms = ^J + J{S-I)-\ (20) 

It is straightforward to check that M5 always is a symmetric matrix: M5 = 
MJ (it suffices for this to use the fact that S^JS = SJS^ = J). 

The symplectic Cayley transform has in addition the following proper- 
ties, which are interesting by themselves: 

Lemma 2 (i) We have 

{Ms + Ms')-^ = -{S' - I){SS' - iy\S - I) J (21) 
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and the symplectic Cayley transform of the product SS' is (when defined) 
given by the formula 

Mss' =Ms + (S^ - I)-^J{Ms + Ms')-^J{S - (22) 

(a) The symplectic Cayley transform of S and are related by 

Mg-i = -Ms. (23) 

Proof, (i) We begin by noting that (20) implies that 

Ms + Ms' = J{I +{S- /)-! + {S' - ly^) (24) 

hence the identity (21). In fact, writing SS' - I = S{S' - I) + S - I, we 
have 

{s' - i){ss' - iy\s -i) = {s' - i){s{s' -i) + s- i)-\s - 1) 

= {{s - i)-^s{s' - i){s' - 1)-^ + {s' - 1)-^)-^ 
= {{s-iy^s+is' -i)-^) 

= I + {S - I)-^ + {S' - I)-^; 

the equahty (21) follows in view of (24). Let us prove (22); equivalently 

Ms + M = Mss' (25) 

where M is the matrix defined by 

M={S'^ - I)-^J{Ms + Ms')-^J{S - ly^ 

that is, in view of (21), 

M = {S^ - I)-^J{S' - I){SS' - 

Using the obvious relations S^ = -JS'^J and (-S"^ +7)-^ = S{S - I)'^ 
we have 

M = {S^ - iy^J{S' - I){SS' - 7)-^ 

= -j{-s-^ + iy\s' - i){ss' - 7)-i 
= -js{s - iy\s' - i){ss' - ly^ 

that is, writing S = S — I + I, 

M = -j{s' - i){ss' - 7)-^ - j{s - iy\s' - i){ss' - iy\ 
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Replacing Ms by its value (20) we have 

Ms + M = 

J{\I +{S- - {S' - I){SS' - I)-^ - (5 - I)-^{S' - I){SS' - 
noting that 

{S - I)-^ - (5 - I)-^{S' - I){SS' - I)-^ = 

{S - i)-\ss' -I-S' + I){SS' - 1)-^) 

that is 

{S - I)-^ -{S- I)-\S' - I){SS' - I)-^ = {S- I)-\SS' - S'){SS' - i)-^ 

= S'{SS' - 1)-^) 

we get 

Ms + M = J{\I - (S" - I){SS' - I)-^ + S'{SS' - I)-^) 
= J(1/ + (5S' -/)-!) 

= Mss' 

which we set out to prove, (ii) Formula (23) follows from the sequence of 
equalities 

Ms-i = \ J + J{S-^ - 
= \ J - JS{S - 
= \J -J{S-I + I){S-I)-^ 
= -Ij - J(S - 1)-' 
= -Ms. 

m 

3.3 Definition and properties of i^(S'oo) 
We define on R^" © R^" a symplcctic form cr® by 

(^^{zi, Z2; z[, 4) = c^{zi,z[) - a{z2, 4) 

and denote by Sp®(2n) and Lag®(2n) the corresponding symplectic group 
and Lagrangian Grassmannian. Let /i® be the ALM index on Lag®, (2n) 
and /xf the Maslov index on Sp® (2n) relative to L G Lag®(2n). 
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For Soo € Sp(^(n) we define 

I.(5oo) = ^/X®((/eS)ooAoo,Aoo) (26) 

where (/ S) oo is the homotopy class in Sp^(2n) of the path 

1 1 — ^ {{z, Stz) : z e M^"} , 0<t<l 

and A = {{z,z) : z eR^""} the diagonal of R^n ^2n_ getting Sf = I® St 
we have Sf G Sp®(2n) hence formulae (26) is equivalent to 

^(^oo) = IfiliSg) (27) 

where /j,^ is the relative Maslov index on Sp®, (2n) corresponding to the 
choice A E Lag®(2n). 

Note that replacing n by 2n in the congruence (6) we have 

f^^iil © 5)00 Aoo, Aoo) = dim((/ 5) A n A) mod 2 

= dimKer(S' - /) mod 2 

and hence 

u{Soo) = ^ dim Ker(5 - /) mod 1. 

Since the eigenvalue 1 of S* has even multiplicity i^(S'oo) is thus always an 
integer. 

The index v has the following three important properties; the third is 
essential for the calculation of the index of repeated periodic orbits (it clearly 
shows that v is not in general additive): 

Proposition 3 (i) For all G Spoo(n) we have 

u{S-') = -u{S^) , u{I^) = (28) 

(loo the identity of the group Sp^{n)). (ii) For all r € we have 

i/(a"5oo) = z^(5oo) + 2r , u{a'^) = 2r (29) 

(Hi) Let Soo be the homotopy class of a path S in Sp(n) joining the identity 
to S E Sp*(n), and let S' £ Sp(n) be in the same connected component 
Sp^(n) as S. Then 1^(5*^) = 1^(5*00) where S'^^ is the homotopy class in 
Sp(n) of the concatenation of and a path joining S to S' in Spo(n). 
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Proof, (i) Formulae (28) immediately follows from the equality {S^)~^ = 
(/©S'~^)oo and the antisymmetry of fi^. (ii) The second formula (29) follows 
from the first using (28). To prove the first formula (29) it suffices to observe 
that to the generator a of 7ri[Sp(n)] corresponds the generator /qq © a of 
7ri[Sp®(2n)]; in view of property (13) of the Maslov index it follows that 

= u{Soo) + 2r. 

(iii) Assume in fact that S and S' belong to, say, Sp'^(n). Let 5*00 be the 
homotopy class of the path S, and S' a path joining S to S' in Sp^(n) (we 
parametrize both paths by t G [0, 1]). Let E^, be the restriction of S' to the 
interval [0,t'], t' < t and 6*00 (^0 the homotopy class of the concatenation 
S * SJ, . We have det(S'(t) - /) > for all t £ [0, f'] hence Sg{t)A n A 7^ as 
t varies from to 1. It follows from the fact that the /i^ is locally constant 
on {Sg : SgA n A = 0} (see §2.3) that the function t 1 — > I^AiSg(.t)) is 
constant, and hence 

^^lisg) = ulisgm = ulisgii)) = ulisg) 

which was to be proven. ■ 

The following consequence of the result above shows that the indices u 
and icz coincide on their common domain of definition: 

Corollary 4 The restriction of the index u to Sp*(n) is the Conley-Zehnder 
index: 

K-Soo) = icz(5oo) if det{S-I)^0. 

Proof. The restriction of u to Sp*(77,) satisfies the properties (CZi), (CZ2), 
and (CZ3) of the Conley-Zehnder index hsted in §3.1; we showed that these 
properties uniquely characterize icz- ■ 

Let us prove a formula for the index of the product of two paths: 

Proposition 5 // ^oo, 5"^, and SooSg are such that det{S — /) / 0, 
det(S' - /) 7^ 0, and det {SS' - /) 7^ then 

u{S^Sg) = u{S^) + u{Sg) + \ sign{Ms + Mg^) (30) 
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where Ms is the symplectic Cayley transform of S; in particular 

u{Sl^) = ru{S^) + i(r - 1) sign Ms (31) 

for every integer r. 

Proof. In view of (27) and the product property (12) of the Maslov index 
we have 

where = I ® S, S"® = 7 © S" and r® is the signature on the symplectic 
space (R2" @ R'^^. a^). The condition det(S'S" - /) 7^ is equivalent to 
S'^S'^A n A = hence wc can apply property (i) in Lemma 1 with i = 
S^S'^A, f = S^A, and £" = A. The projection operator onto S^S'^A 
along A is easily seen to be 



Prgeg/eA.A 



{i-ss')-^ -{i-ss')-^ 
SS'{I - ss')-^ -SS'{I - ss')-^ 



hence r®(S'®S"® A, 5® A, A) is the signature of the quadratic form 

Q{z) = (7©(Pr5e5'eA,A(2, Sz); {z, Sz)) 

that is, since o"® = a Q a: 

Q{z) = a{{I - SS')-\I - S)z, z)) - a{SS'iI - SS')-\l - S)z, Sz)) 
= a{{I - SS')-Hl - S)z, z)) - a{S'{I - SS')-\l - S)z, z)) 
= ai{I-S'){I-SS')-\l-S)z,z)). 

In view of formula (21) in Lemma 2 we have 

(7 - S'){SS' - I)-\I -S) = {Ms + Ms')-^J 

hence 

Q{z) = -{{Ms + Ms')-'Jz,Jz) 
and the signature of Q is thus the same as that of 

Q'{z) = -{{Ms + Ms')-h,z) 

that is — sign(Ms + Mg/). This proves formula (30). Formula (31) follows 
from (30) by induction on r. ■ 
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It is often deplored in the literature on Gutzwiller's formula (1) that it is 
not always obvious that the index of the periodic orbit 7 is independent 
on the choice of the origin of the orbit. Let us prove that this property 
always holds: 

Proposition 6 Let (ft) be the flow determined by a (time-independent) 
Hamiltonian function on M^*^ and z such that friz) = z for some 
T > 0. Let z' = ft'{z) for some t' and denote by St{z) = Dfxiz) and 
St{z') = Dfxiz') the corresponding monodromy matrices. Let St{z)oo o,nd 
St{z')oo be the homotopy classes of the paths t 1 — > St{z) = Dft{z) and 
t ^ St{z') = Dft{z'), 0<t<T. We have u{St{z)^) = u{St{z'U). 

Proof. The monodromy matrices St{z) and St{z') are conjugate of each 
other; in fact (proof of Theorem 6 in [10]): 

ST{z') = St'{z')ST{z)St'{z'r'; 

since we will need to let t' vary we write St{z') = St{z' ,t') so that 

St{z' ,t') = St,{z')ST{z)St>{z')-\ 

The paths 1 1 — > St{z') and 1 1 — > St'{z')St{z)St'{z')'~^ being homotopic with 
fixed cndpoints St{z' ,t')oo is also the homotopy class of the path t 1 — > 
St'{z')St{z)St'{z')-^. We thus have, by definition (26) of 

v{ST{z\t')^) = \lil^XS^{zU) 

where we have set 

Aj/ = (/ e St'{z')-^)/^ and S^{z)oo = L^® St{z)^. 

Consider now the mapping t' 1 — > ji^^^{Sj-{z)(xi)\ we have 

S^{z)At>f\l^t' = {z-Sz = z} 

hence the dimension of the intersection S'^{z)^ti fl A^/ remains constant as 
t' varies; in view of the topological property of the relative Maslov index the 
mapping t' 1 — ^ iJ,'^^^{S^{z) 00) is thus constant and hence 

I^{STiz',t'U) = 1^{St{z',0)^) = iy{ST{z)oo) 

which concludes the proof. ■ 
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3.4 Relation between v and /x^^ 



The index v can be expressed in simple - and useful - way in terms of the 
Maslov index jj^p on Sp(^(n). The following technical result will be helpful 

in establishing this relation. Recall that S G Sp(n) is said to be "free" if 
we have Sip H £p = 0; this condition is equivalent to det B ^ Q when S is 
identified with the matrix 

\A B 
C D 



(32) 



in the canonical basis. The set of all free symplectic matrices is dense in 
Sp(n). The quadratic form W onW^ x defined by 

W{x, x') = ^ {Px, x) - {Lx, x') + ^ {Qx', x') 

where 

P = DB-\ L = B-\ Q = B-^A (33) 
then generates S in the sense that 

{x,p) = S{x',p') ^p = d^W{x,x') , p' = -d^>W{x,x'). 

We have: 

Lemma 7 Let S = Sw G Sp(ra) be given by (32). We have 

det{Sw -I) = {-^T det B det{B-^A + DB'^ - B'^ - (S^)"^) (34) 
that is: 

det{Sw-I) = (-l)"det(L-^)det(P + Q-L-L^). 
In particular the symmetric matrix 

P + Q-L-L^ = DB-^ + B-^A - B'^ - (S^)"^ 
is invertible. 

Proof. Since B is invertible we can factorize S — I as 



A- 1 B 
C D-I 



B 

1 D-I 



'C-{D-I)B-\A-I) 
B-\A-I) I 
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and hence 



det{Sw -I) = dct(-B) dct(C - {D - I)B-\A - I)) 

= {~lYdctBdct{C -{D-I)B-^{A-I)). 

Since S is symplectic we have C — DB~^A = —{B^)~^ and hence 

C-{D- I)B-^{A - I)) = B-^A + DB-^ - B'^ - (5^)"^ 

the Lemma follows. ■ 

Let us now introduce the notion of index of concavity of a Hamiltonian 
periodic orbit 7, defined for < t < T, with 7(0) = ^yiT) = Zq. As t goes 
from to T the linearized part D^{t) = St{zQ) goes from the identity to 
St{zq) (the monodromy matrix) in Sp(n). We assume that St{zo) is free 
and that det(S'T(2;o) - /) / 0. Writing 



A{t) Bit) 
Cit) Bit) 



we thus have dei B{t) 7^ in a neighborhood [T — £,T + e] of the time T. 
The generating function 

W{x,x',t) = ^ {P{t)x,x) - {L{t)x,x') + ^ {Q{t)x',x') 

(with P{t), Q{t), L{t) defined by (33) thus exists for T - e < t < T + £. By 
definition Morse's index of concavity [17] (also see [18, 19]) of the periodic 
orbit 7 is the index of inertia 

Inert W^^ = Inert(P + Q - L - L^) 

of W^'a;, the matrix of second derivatives of the function x 1 — > W{x,x;T) 
(we have set P = P{T), Q = Q{T), L = L{T)). 

Let us now prove the following essential result; recall that denotes 
the reduced Maslov index (15) associated to jif. 

Proposition 8 Let t 1 — > St be a symplectic path, < t < 1. Let Soo & 
SPoo(^) be the homotopy class of that path and set S = Si. If det{S — I) ^ 
and Sip n £p = then 

J^iSoo) = ^il^epiSoo) + sign W^J = mep{Soo) - Inert W^^ (35) 

where Inert W^^ is the index of concavity corresponding to the endpoint S of 
the path t 1 — > St- 
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Proof. We will divide the proof in three steps. Step 1. Let L G Lag®(4n). 
Using successively formulae (27) and (16) we have 

^^(^oo) = liii2iSg)+r^{S'^A,A,L)-T^iSQA,S^L,L)). (36) 
Choosing in particular L = Lq = £p (B £p we get 

/.e (5© ) = //©((/ e suiip e £p), {£p e ip)) 

= /^(^P,oo, ^P,oo) — l^{ip,oo, Sooip,oo) 
= -IJ'{^P,oo,Sooip,oo) 

= Wp('S'oo) 

so that there remains to prove that 

T^{S^A,A,Lo)-r^{S^A,S^Lo,Lo) = -2sign<,. 

Step 2. We are going to show that 

re(5©A,,SeLo,Lo) = 0; 

in view of the symplectic invariance and the antisymmetry of r® this is 
equivalent to 

r©(Lo,A,Lo,(Se)-%) = 0. (37) 

We have 

AnLo = {(0,p;0,p) :pGM"} 

and (5©)-^LonLo consists of all (0,^,5-^(0,/)) with ^-^(O,^') = (0,^); 
since S (and hence S~^) is free we must have p' = p" = so that 

(5©)-iLo n Lo = {{0,p; 0,0) -.pe W}. 

It follows that we have 

Lo = A n Lo + {S'^y^Lo n Lo 

hence (37) in view of property (ii) in Lemma 1. Step 3. Let us finally prove 
that. 

re(5eA,A,Lo) = -2signW^^;; 

this will complete the proof of the proposition. The condition det(6' — /) 7^ 
is equivalent to S^A H A = hence, using property (i) in Lemma 1: 

T^{S^A,A,Lo) = -T^{S^A,Lo,A) 



20 



is the signature of the quadratic form Q on Lq defined by 

Q(0,p,0,p') = -c7©(Pr5eA,A(0,P,0,p');0,P,0,p') 

where 

^^""^'^ s{s-i)-^ -~s{s^ir^ 



is the projection on S^A. along A in M "-©M It follows that the quadratic 
form Q is given by 

Q(0,p,0,p') = -^®((^-5)"'(0,/),.S(/-5)-^(0,/);0,p,0,p') 

where we have set p" = p — p'; by definition of cr® this is 

Q(0,p,0,y) = -a{{I - S)-\Qy ), (0,p)) + a{S{I - S)-Ho,p"), {0,p')). 

Let now Mg be the symplectic Cayley transform (19) of S; we have 

{I - S)-^ = JMs + p , S{I-S)-^ = JMs-\l 

and hence 

Q(0,p,0,i>') = -a{{JMs + \my),{Q,p)) + a{{JMs - \l){Qy),M) 
= -a{JMs{0,p"), (0,p)) + a{JMs{0,p"), {0,p')) 
= a{JMs{0,p"),{0,p")) 
= -{Ms{0,p"),{0,p")). 



Let us calculate explicitly M5. Writing S in usual block- form we have 
S-I 



B 

1 D-I 



C - {D - I)B-\A- I) 
B-\A-I) I 



that is 



B 

1 D-I 



W" 







B-^{A-I) I 



S-I 

where we have used the identity 

C-{D- I)B-\A - I)) = B-^A + DB-^ - B'^ - {B^^Y^ 
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which follows from the relation C — DB = —{B^) ^ (the latter is a 
rephrasing of the equalities A — B^C = I and B = B^D, which 
follow from the fact that S^JS = S^JS since S G Sp(n)). It follows that 



{s-iy 







(I - D)B-^ I 
B-^ 



Ms 



B~\I - A)(W^,)-^ I 

{W'^x)-\I - D)B-^ «J-i 
B-\I - A){W'^,)-\I - D)B-' + B-^ B-\I - A)KJ-i 

and hence 

B-\I - A){W^^)-\I - D)B-' + B-' \I + B-\I - A){W'^,)-^ 
-ll - {W^^)-\I - D)B-^ -W.r 

from which follows that 

Q{Q,p^Q,p') = {{w':,)-v y) 

= {{w':,)-\p-p\{p-p')). 

The matrix of the quadratic form Q is thus 



and this matrix has signature 2 sign(VF"^)^^ = 2 sign H^"^, proving the first 
equality (35); the second equality follows because //^p(S'oo) = 2m£p(>S'oo) — n 
since Sip Dip 



and the fact that W"^ has rank n in view of Lemma 7. 



Remark 9 Lemma 7 above shows that if S is free then we have 
— arg det(S' — I) = n + arg dct B + arg det W^^ mod 2 

TT 

= n — arg det B + arg det W"^ mod 2 

In [5, 6] we have shown that the reduced Maslov index m£p(Soo) corresponds 
to a choice of arg det B modulo 4; Proposition 8 thus justifies the following 
definition of the argument of det{S — I): 

— argdet(5 — I) = n — I'iSoo) mod 4. 

TT 

That this is indeed the correct choice m,odulo 4 has been proven by other 
means ( the Weyl theory of the metaplectic group ) by one of us in a recent 
publication [11]. 



22 



4 An Example 



Let us begin with a very simple situation. Consider the one-dimensional 
harmonic oscillator with Hamiltonian function 



H 



all the orbits are periodic with period 27r/a;. The monodromy matrix is 
simply the identity: = / where 



cos ut sin ijjt 
— sin ut cos ut 



Let us calculate the corresponding index i/(Eoo). The homotopy class of path 
1 1 — > St as t goes from to T = 27r/a; is just the inverse of a, the generator 
of 7ri[Sp(l)] hence z/(Sqq) — — 2 in view of (29). If we had considered t 
repetitions of the orbit we would likewise have obtained viTioo) = — 2r. 

Consider next a two-dimensional harmonic oscillator with Hamiltonian 
function 



H 



we assume that the frequencies Uy, Ux are incommensurate, so that the only 
periodic orbits are librations along the x and y axes. Let us focus on the orbit 
7a; along the x axis; its prime period is T = 21^ /ux and the corresponding 
monodromy matrix is 



Si 



1 0' 

cosx sinx 

10 

—sinx cosx 



it is the endpoint of the symplectic path t 
the matrices 



St, < t < 1, consisting of 



cos27rt sin27rt 

cosxi sinxi 

— sin 2Trt cos 27rt 

— sin xt cos xt 

In Gutzwiller's formula (1) the sum is taken over periodic orbits, including 
their repetitions; we arc thus led to calculate the Conley-Zehnder index of 
the path 1 1 — > St with < t < r where the integer r indicates the number of 
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repetitions of the orbit. Let us calculate the Conley-Zehnder index z^(S'r,oo) 
of this path. We have S't = St where 



cos 2Trt sin 27rt 
— sin 27rt cos 27rt 



, i^t 



cos xt sin xt 
— sin xt cos xt 



in view of the additivity property of the relative Maslov index we thus have 
where the first term is just 

in view of the calculation wc made in the one-dimensional case with a dif- 
ferent parametrization. Let us next calculate i^{Sr,cx>)- We will use formula 
(35) relating the index u to the Maslov index via the index of concavity, so 
we begin by calculating the relative Maslov index 

mip{Sr,oo) = ?^('5r,oo^P,oo,^Poo)- 

Here is a direct argument; in more complicated cases the formulas we proved 
in [10] are useful. When t goes from to r the line Stip describes a loop 
in Lag(l) going from £p to Srip- We have St G U{1); its image in C/(1,C) 
is e~*^* hence the Souriau mapping identifies Stip with e~^*^*. It follows, 
using formula (9), that 

me,{Sr,^) = ^ {-2rx + ^ Log(-e-2-x)) + ^ 

= ^(-2rx + iLog(e*(-2'->^+-)))+^ 

The logarithm is calculated as follows: for 9 ^ (2k + l)7r {k € Z) 

+ 7r" 



Loge'" = ie -2Tri 



27r 



and hence 



Log(e^(-2''^+'')) = -i(2rx + tt + 27r 
it follows that the Maslov index is 

'rx 



rx 



L TT 



mip{Sr,oo) 



TT 



(38) 
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Inert W"^ 



To obtain i'{Sr,oo) we note that by (35) 

'^('S'r,oo) = "T-fp (5*1,0 

where Inert W^^ is the concavity index corresponding to the generating func- 
tion of Sf, the latter is 

1 



W{x,x',t) 



{{x^ + x'^) cos xt - 2xx') 



2 sin xt 

hence W"^ = — tan(xt/2). We thus have, taking (38) into account, 

'rx] 



Inert 



, rx 
tan — 
2 



a straightforward induction on r shows that this can be rewritten more 
conveniently as 

l^{Sr,oo) = -1 - 2 

Summarizing, we have 

[rx] 



rx 
27r 



-2r 



2tt 



rx 



hence the index in Gutzwiller's formula corresponding to the r-th repetition 
is 

IJ'X,r = -l'{Sr,oo) = 1 + 2r + 2 

that is, by definition of x, 



2-K 



1 + 2r + 2 



confirming the calculations in [1, 4, 20, 23]. 

Remark 10 The calculations above are valid when the frequencies are in- 
commensurate. If, say, ujx = uJy, the calculations are much simpler: in this 
case the homotopy class of the loop t i — > 5"^, < i < 1, is © and 
by the second formula (29), 



which is zero modulo 4 (cf. [20]). 



4r 



Acknowledgement 11 One of the authors (M. de Gosson) has been sup- 
ported by the grant 2005/51766-7 of the Fundagdo de Amparo d Pesquisa do 
Estado de Sao Paulo (FAPESP); it is a pleasure for him to thank P Piccione 
for his kind hospitality. 



25 



References 

M Brack and R K Bhaduri 1997 Semiclassical Physics (Addison- 
Wesley) 

S E Cappell, R Lee and E Y Miller 1994 Comm. Pure and Appl. Math. 
17 121 

C Conley and E Zelinder 1984 Comm. Pure and Appl. Math. 37 207 

S C Creagh, J M Robbins, and R G Littlejohn 1990 Phys. Rev. A, 42(4) 
1907 

M de Gosson 1990 Ann. Inst. Fourier 40(3) 537 

M de Gosson 1992 J. Geom. Phys. 9 255 

M de Gosson 1992 J. Math. Pures et Appl. 71 429 

M de Gosson 1997 Maslov Classes, Metaplectic Representation and La- 
grangian Quantization {Research Notes in Mathematics 95) (Berlin: 
Wiley-VCH) 

M de Gosson 2001 The Principles of Newtonian and Quantum Me- 
chanics: the Need for Planck 's Constant h; with a foreword by B. Hiley 
(London: Imperial College Press) 

M de Gosson and S de Gosson 2003 J. Phys A: Math, and Gen. 36 615 

M de Gosson 2005 Letters in Mathematical Physics 72 129 

M C Gutzwiller 1990 Chaos in Classical and Quantum Mechanics (In- 
terdisciplinary Applied Mathematics) (Berlin: Springer) 

H. Hofer, K Wysocki, and E Zehnder 1995 Geometric and Functional 
Analysis 2(5) 270 

J Leray 1981 Lagrangian Analysis and Quantum Mechanics, a mathe- 
matical structure related to asymptotic expansions and the Maslov index 
(Cambridge, MA: MIT Press, ) 

G Lion, and M Vergne 1980 The Weil representation, Maslov index and 
Theta series {Progress in mathematics 6) (Boston: Birkhauser) 

B Mehlig and M Wilkinson 2001 Ann. Phys. 18(10), 541 



26 



[17] M Morse 1935 The Calculus of Variations in the Large (Providence, R 
I: American Mathemtical Society) 

[18] P Muratore-Ginnaneschi 2003 Phys. Rep. 383 299 

[19] R C Nostre Marques, P Piccione, and D Tausk 2001 in Differential 
Geometry and its Applications (Opava: Conference Proceedings) 

[20] M Plctyukhov and M Brack 2003 J. Phys A: Math, and Gen. 36 9449 

[21] J M Robbin 1991 Nonlinearity 4 343 

[22] J M Souriau 1975 In: Group Theoretical Methods in Physics, Lecture 
Notes in Physics, Springer- Verlag, 50 17 

[23] A Sugita 2001 Ann. Phys. 288 277 

[24] C T C Wall 1969 Invent. Math. 7 269. 



27 



